A road to hydrogenating graphene by a reactive ion etching plasma 
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We report the hydrogenation of single and bilayer graphene by an argon-hydrogen plasma pro- 
duced in a reactive ion etching (RIE) system. Electronic transport measurements in combination 
with Raman spectroscopy are used to link the electric mean free path to the optically extracted 
defect concentration. We emphasize the role of the self-bias of the graphene in suppressing the 
erosion of the flakes during plasma processing. We show that under the chosen plasma conditions 
the process does not introduce considerable damage to the graphene sheet and that hydrogenation 
occurs primarily due to the hydrogen ions from the plasma and not due to fragmentation of water 
adsorbates on the graphene surface by highly accelerated plasma electrons. For this reason the 
hydrogenation level can be precisely controlled. The hydrogenation process presented here can be 
easily implemented in any RIE plasma system. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 

Hydrogenation of carbon materials, e.g. graphite, car- 
bon nanotubes or carbon foams, has triggered a large 
technological and scientific interest, with its main focus 
on hydrogen physisorbtion in hydrogen storage systems 
However, for electronic applications the chemisorp- 
tion of hydrogen is even more interesting, as it allows 
for tuning of electronic properties in carbon conjugated 
systems. An excellent candidate for such manipulation 
is graphene, a single layer of graphite, built only from 
sp^ carbons and demonstrating high carrier mobilities 
0. Similarly to single wall nanotubes Q, the small vol- 
ume and large contact area of graphene makes chemisorp- 
tion of hydrogen an efficient way to modify its electronic 
properties [1,1^. Depending on the H coverage one can 
tune the transport properties of graphene from metallic 
to semiconducting, and ultimately to an insulating state 
for its fully hydrogenated derivative graphane i6] . Open- 
ing of a bandgap by hydrogenation in otherwise gapless 
graphene can be also an elegant way to fabricate a cir- 
cuit consisting of a single material: graphene, with both 
metallic and semiconducting parts. 

Apart from microelectronic applications, the influence 
of hydrogen on electronic transport in graphene has great 
scientific relevance as well. In particular to understand 
the role of localized defects as scattering centers limit- 
ing carrier mobility 0], the transition in charge trans- 
port from the Drude type (in pristine graphene) to the 
variable range hopping type (in strongly hydrogenated 
graphene) [6| , or in predictions of magnetism originating 
from hydrogen defects [S, 9] . 

Chemisorption of hydrogen on a graphene surface 
changes the carbon electronic orbitals from sp^ to sp'^ 
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hybridization and results in a localized state. The po- 
tential barrier for hydrogen adsorption to the surface of 
graphene is about 0.2 eV [1, Part of this energy 

is consumed by the displacement of the carbon out from 
the graphene plane to obtain the tetragonal sp^ geometry. 
This adsorbtion barrier is lower in initially curved or pro- 
truding structures (at grain boundaries, lattice defects 
or on ripples), where structural deformation is already 
present For effective and controllable hydrogena- 

tion of graphene several techniques have been explored 
so f ar, including exposure to an atomic hydrogen source 
pT| - [l3| . electron beam (e-bearn) exposure of highly hy- 
drated lithography resist HSQ [iJl and e-beam exposure 
of a water adhesive layer on graphene [l^ [l^ . Among 
this techniques exposure to an argon-hydrogen plasma 
produced in a DC ^ or RF source seems promising 
alternative due to the high energy and reactivity of the 
incident hydrogen ions, what enables their chemisorption 
even on the flat surface of graphene. Hydrogenation by 
an Ar/H2 plasma can lead to a high and fast hydrogen 
uptake, it does not require special sample preparation 
and is compatible with microfabrication techniques. 

Estimation of the H content in micromechanically 
cleaved graphene flakes after hydrogen treatment is very 
difficult. The standard methods known from graphite, 
like thermal programmed desorption (TPD) _18;-20j , are 
insensitive to the possible amounts of desorbed hydrogen 
from micro-sized fiake. Estimation of H coverage from 
scanning tunneling microscopy (STM) topography im- 
ages carries limitations that STM probes the surface only 
locally, measurements are time consuming and difficult 
when graphene is deposited on the insulating substrate. 
An appealing alternative is Raman spectroscopy, which 
is a relatively easy, non-destructive, non-contacting and 
quick method to probe H coverage from even microme- 
ter sized samples and can be carried out at room tem- 
perature and atmospheric pressure. Chemisorption of H 
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induces Raman bands which are normally symmetry for- 
bidden in the graphene spectrum. The assignment of 
these bands to hydrogen adsorbates allows an indirect 
estimate of the H content [2l| . 

In this work we demonstrate the hydrogenation of 
graphene by an RF plasma of an argon-hydrogen gas 
mixture using reactive ion etching (RIE). This technique 
has not been explored for graphene hydrogenation so far, 
despite the fact that RIE is widely used for electronic 
device microfabrication. We characterize the hydrogena- 
tion properties of the RF plasma and its reversibility 
under moderate thermal annealing by means of Raman 
spectroscopy. Further we present the electronic transport 
measurements in single layer (SLG) and bilayer graphene 
(BLG) which enables us to relate the structural defects to 
graphene transport properties. In the control experiment 
we compare the effect of the Ar/H2 plasma with the pure 
Ar plasma in two types of sample (in bare flakes on insu- 
lating substrate and in graphene devices). The observed 
differences highlight the role of the floating potential of 
the non-contacted graphene flakes for acceleration of the 
graphene erosion. As this effect is completely suppressed 
in graphene devices, we conclude that there graphene hy- 
drogenation happens primarily due to the hydrogen ions 
and not to highly accelerated plasma electrons fragment- 
ing water add-layer on graphene surface, as suggested in 
Ref. [22|] . 
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FIG. 1: (Color online) Schematic drawing of the reactive ion 
etching chamber used in the experiment. The sample 
is placed on the bottom electrode, connected to the 
RF power source, while the top electrode and the 
walls are grounded. On the right the corresponding 
potential profile during the plasma operation. The 
positive bias Vp inside the plasma and negative self- 
bias N SB at the sample are indicated. Graphene can 
be insulated from the bottom electrode by the Si02 
substrate or connected to it by fabricated metallic 
contacts. In latter case graphene will have the same 
self-bias as the bottom electrode. 



II. EXPERIMENTAL METHODS 
A. RF plasma conditions 

The hydrogenation is performed in a reactive ion etch- 
ing reactor with a parallel plate geometry, schematically 
depicted in Fig. [T] The diameter of the bottom elec- 
trodes, on which the samples are placed, is 300 mm while 
the opposite top wall of the chamber serves as a grounded 
counter-electrode. A high frequency generator operat- 
ing at 13.56 MHz is capacitively coupled to the bottom 
electrode, and a matching of the electrical network to 
the plasma is accomplished by mechanical tuning of the 
impedance in the circuit. 

In view of safety considerations, we use a gas mixture 
of H2 (15%) with Ar (85%) as a balance gas. The ioniza- 
tion energy of Ar, YjAt = 15.76 eV, is very close to the 
ionization energy of II2, E//^ = 15.42 eV, therefore the 
induced plasma is composed of ions from both species. 
The inlet of gas is controlled by an Ar mass flow con- 
troller. In all presented plasma hydrogenation processes 
the gas flow is kept constant at 200 seem and the pressure 
in the chamber is 0.05 mbar. To reduce the reactivity of 
the plasma, especially carbon sputtering by Ar ions, we 
use the plasma at the lowest ignition power, P = 3 W 
(power density is ~ 4 mW/cm^), and we tune the circuit 
impedance to reduce the built-in DC self-bias between 
the bottom electrode and the plasma (Vss in Fig. [1]) 
down to zero. Reducing Ysb to zero minimizes the ion 



acceleration and possible sputtering effects on graphene. 
We analyze two cases, one where the graphene flakes are 
electrically insulated from the chamber electrodes (by the 
Si02 substrate) and one where the flake is in electric con- 
tact with the source electrode (latter called a graphene 
device). In the first case, the potential of the flake is float- 
ing, which may result in negative charging of the flake 
before the plasma quasi-equilibrium state and Vsb=0 
is reached (in the first 3 seconds after the plasma igni- 
tion). This charging is largely suppressed for graphene 
device, which is in electrical contact with the chamber 
electrode. On the basis of work of Nunomura et al. 
we estimate that we are in a collisional regime, with ion 
bombardment energy in the range of 5-20 eV and that 
the dominant hydrogen radicals are , with a much 
smaller concentrations of H2 and . We note that the 
processing conditions and hydrogenation speed are diffe- 
rent from the one explored in Ref. [l^ ■ The gas pressure 
in that process is 2 orders of magnitude higher than it is 
here, and Luo et al. used a grounded bottom electrode. 



B. Raman spectroscopy of prisitine and 
hydrogenated graphene 

Information about the H content can be obtained in- 
directly from Raman spectra [iTj . In pristine graphene 
only two vibrational modes are Raman active: an in- 
plane optical vibration of E2g symmetry, at 1580 cm"-'^. 
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called G band, produced by sp^ carbon network and 
a resonantly enhanced two phonon scattering process, 
around 2670 cm~^, called 2D (or sometimes G'). The 
presence of sp^ defects breaks the translational symme- 
try in the graphene lattice and activates other resonant 
transitions. The most significant is so-called defect band 
D at 1340 cm~^, forbidden in the ideal sp^ graphene lat- 
tice. The D band results from a second order process in- 
volving intervalley elastic scattering of electron by defect 
and inelastic scattering on phonon. It is worth noting 
that the 2D mode is an overtone of the D peak, with 
the difference that in case of the 2D band an electron is 
scattered by a second phonon instead of a defect. Ad- 
ditionally sp^ defects induce a much weaker D' band at 
1620 cm~^, coming from intravalley defect scattering and 
a peak, which can be assigned to the combination of the 
D and G mode (G+D) at - 2940 cm'^ [2l|. These prop- 
erties of graphene make Raman spectroscopy a sensitive 
tool for detection of chemisorbed H defects. It is worth 
noting that the physisorbed molecules do not change the 
hybridization of carbons and hence do not contribute to 
the Raman signal of the D band. In the Ar/H2 plasma 
process presented here, one has to take into account also 
the effect of the Ar ions, which by bombarding graphene 
could induce other sp'^-type defects: vacancies. These 
defects also contribute to the D band intensity in Raman 
spectra; therefore, care must be taken when one assigns 
the D band intensity solely to the H adatoms. Later in 
this work we prove by studies of thermal desorption that 
the sputtering effect of Ar ions is largely suppressed in the 
chosen plasma conditions (considerably low RF power, 
high gas pressure). This assures the hydrogen origin of 
sp^ defects. To quantify the level of hydrogenation we 
use the integrated intensity ratio / Iq of Raman bands, 
which relates the amount of sp'^ defects in the graphene 
lattice to its inherent sp^ bonds. Raman spectra are ob- 
tained using a Horiba T64000 micro-Raman spectrom- 
eter with 532 nm laser excitation wavelength, spectral 
resolution of ^ 2 cm~^, laser spot size <10 /im in di- 
ameter and power density below 0.5 mW to avoid laser 
induced heating. First, we study the evolution of the D 
band and its amplitude in comparison with the G band 
in Raman spectra at various plasma exposure times. For 
that purpose we select a set of graphene flakes deposited 
on Si02/Si substrate (300 nm of Si02) by micromechan- 
ical cleavage of Kish graphite. For each flake we obtain 
a pristine Raman spectrum. With that we exclude the 
presence of initial disorder. By analyzing the shape and 
FWHM of the 2D band we confirm the number of layers 
in the chosen flakes [l^, . Then each sample is exposed 
separately to the Ar/Il2 plasma for a specific amount of 
time and immediately after that the Raman spectrum in 
ambient conditions is acquired. 

A typical Raman spectrum of graphene before and af- 
ter the plasma exposure is shown in Fig. [2al Hydrogena- 
tion results in activation of additional vibrational modes, 
two of which are depicted in Fig. I2al a disorder induced 
D band at ~ 1340 cm"! and D' band at 1620 cm'^ The 
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FIG. 2: (Color online) (a) Raman spectrum of pristine single 
layer graphene (black) and after 20 min of exposure 
to the Ar/H2 plasma (blue). Exposure induces ad- 
ditional Raman bands: a D band around 1340 cm~^ 
and a weaker D' band around 1620 cm~^. The in- 
crease of FWHM of original graphene bands (G, 2D) 
is apparent, (b) Integrated intensity ratio between 
the D and G bands of SLG after different Ar/H2 
plasma exposure times. The scattering of the data 
for different samples is attributed to the floating po- 
tential of the graphene flake during exposure, (c) 
The change of the Id/Ig ratio of exposed flakes un- 
der annealing on hot-plate for 1 min. The plasma 
exposure time for each flake is indicated next to the 
corresponding Id/Ig values. In flakes exposed for 
less than 1 h the D band could be almost fully sup- 
pressed {Id/Ig <0.2), which confirms the H-type ori- 
gin of defects. In longer-exposed samples (80 min and 
2 h) annealing does not significantly reduce Id/Ig, 
which suggests a different nature of defects there, e.g. 
vacancies. 



evolution of the integrated intensity ratio between the D 
and G bands in the Raman spectrum, /Iq , after diffe- 
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rent exposure times is shown in Fig. l2bl Wc note its sim- 
ilar behavior to that presented in Ref . JTJ . The increase 
of the exposure time results in the increase of the ratio be- 
tween the D and G band up to the point where there are 
so many defects in the graphene lattice that the graphene 
electronic band structure is degraded, reducing possible 
optical transitions for both D and G bands [26j . The ini- 
tial increase and then decrease of the Id / Ig ratio with 
an increasing number of defects in graphene is reported 
irrespective of the origin of the defects |27l - [30| . After hy- 
drogenation, all original Raman bands of graphene show 
an increase of their FWHM, which is attributed to the 
local deformation of the lattice and a larger variation in 
vibrational/phonon energy. 



C. Reversibility of hydrogenation under annealing 

To confirm that the defects in graphene detected by 
Raman spectroscopy originate from H adsorbates, we 
study the change of the Id/Ig ratio after heat treat- 
ment. The comparative studies of hydrogen desorption 
in graphite by TPD show that H starts to desorb already 
at moderate temperatures, >100 °G, with the desorption 
maxima at 175 °G and 290 °G and estimated activation 
energy for desorption is 0.6 eV Note that these 

temperatures are too low to heal possible vacancies in 
graphene. We perform the heating in a nitrogen envi- 
ronment on a hot-plate, with temperature ranging from 
75 °C to 275 °C, each time for 1 min. As can be seen in 
Fig. [2cl heating results in a decrease of the ratio Id/Ig- 
It starts already at 75 °G and continues decreasing with 
increase of heating temperature. Desorption of hydrogen 
below 100 °C, also reported in can originate from 
different nature of hydrogenation by plasma in compar- 
ison with an atomic hydrogen source, as energetic ions 
can bind to graphene in more diverse, also meta-stable, 
configurations of hydrogen clusters. 

After heating at 275 °G, Id/Ig drops below 0.2 in the 
case of the samples exposed to plasma for less than 1 h. 
The samples exposed for 80 min and 2 h show a much 
smaller decrease of defect band intensity with tempera- 
ture. This means that after prolonged exposures the D 
band in these flakes must originate primarily from car- 
bon vacancies rather than H adsorbates. In a control 
experiment we expose the graphene flakes to the pure Ar 
plasma at the same RIE exposure conditions. We observe 
that the pure Ar plasma induces substantial etching of 
graphene, with complete erosion of the flake after about 
30 min. The different etching rate of the Ar/H2 plasma 
versus pure Ar plasma can be explained by the mass dif- 
ference between H and Ar ions. Lighter H ions are faster 
accelerated by the bias difference between the plasma 
and graphene and they reach the graphene surface sooner 
than Ar ions. By charge transfer H ions effectively neu- 
tralize the negative potential of the flake, reducing the 
self-bias voltage between the sample and the plasma and 
the acceleration of much heavier Ar ions. Although the 



carbon vacancies seem to contribute substantially to the 
D band signal after the plasma exposures with Ar, this ef- 
fect is completely suppressed in graphene devices, where 
the flake is in electric contact with bottom electrode. Ex- 
posure of the contacted flake to the Ar plasma did not 
produce any defect related Raman bands even after pro- 
longed exposure (>3 h) and later in the text we show no 
significant change in graphene electronic mobility under 
the Ar plasma exposure. This emphasizes the role of the 
floating potential of the graphene sample for amplifying 
the etching speed. 



D. Electronic transport in hydrogenated graphene 

To gain more information about the role of different H 
coverage on electronic transport, we perform 4 terminal 
resistivity measurements in single and bilayer graphene 
devices after sequential exposure to the Ar/H2 plasma 
(devices are exposed simultaneously) . The measurements 
are done at room temperature and in vacuum shortly af- 
ter the plasma exposure. The inset of Fig. [5a] shows ex- 
emplary resistivity measurements at different charge car- 
rier concentrations for SLG device. The carrier concen- 
tration n can be extracted from the charge induced by the 
gate voltage Vg with respect to the voltage of the charge 
neutrality point (CNP) Vd (also called Dirac point, 
where the valence band of graphene touches the conduc- 
tion band) by using the formula: n — Cg/e{VD — Vg), 
where gate capacitance Cg=115 aF/firn^ for 300 nm 
Si02. Upon exposure the position of the Dirac point 
shifts to positive voltages, indicating the hole doping 
from H. Linking this shift directly to the amount of ad- 
sorbed H is however not appropriate here, as the mea- 
surements are done ex-situ and other dopants, like ph- 
ysisorbed water molecules, could screen the doping in- 
duced by H [Slj . For that reason we focus on the resisti- 
vity changes at the charge neutrality point and in a high 
doping regime, where graphene shows metallic behavior 
(here arbitrarily taken at ^ 2 x lO^'^ cm^'^). In Fig. [3a| 
one can see that with the increase of the exposure time 
the SLG resistivity changes from a few kil to Mil and for 
BLG to hundreds of kfi. Upon hydrogenation the resis- 
tivity difference between CNP and a high doping regime 
changes from ^^3 kfi to ^300 kSl, and its gate voltage 
characteristic broadens indicating the large amount of 
charge impurities/inhomogeneities. (If one defines the 
width of resistivity dependence p from the charge car- 
rier concentration as the distance between its deflection 
points, then upon hydrogenation this width changes in 
SLG from 8 x 10" cni^^ to > 1 x 10^"* cm^^)^ ^s one 
might expect, the increase of graphene resistivity with ex- 
posure time is slower for BLG than for SLG, as there the 
graphene layer underneath is unexposed. Moreover, BLG 
shows a monotonic increase of resistivity with exposure, 
whereas for SLG we observe a non-monotonic change 
in resistivity, which suggests a change in the transport 
mechanism for exposure times >30 min. The same be- 
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havior is reflected by the electron mean free path I, cal- 
culated here using the formula: / = 2D/vp, where vp is 
Fermi velocity of electrons in graphene, vp = 10^ m/s, 
and D is a diffusion coefficient (obtained from Einstein 
relation D = uje^v^ v is the density of states). In 
the calculation we neglect the effect of finite tempera- 
ture on the density of states (DOS) and any broaden- 
ing due to charge impurities; the interlayer coupling in 
DOS of bilayer graphene 71 = 0.4 eV, after [33]). Fig- 
ure I3bl shows a change of the mean free path I with the 
H plasma exposure. It decreases monotonically for BLG 
and non-monotonically for SLG. The shaded area marks 
the mean free path distances below the length of the C- 
C bond (^1.4 A), where the diffusion transport model 
loses its physical meaning. The fact that the estimated 
mean free path for SLG after '^2 h of exposure enters 
this range provides us with evidence that the transport 
there can no longer be described by the semi-classical 
Drude model. Low temperature measurements presented 
in Ref. Q show that in the heavily hydrogenated samples 
the transport enters a variable range hopping regime, but 
the full description of this transition is still lacking. 

Additionally, in a control experiment we perform the 
same electrical characterization of graphene devices ex- 
posed to the pure Ar plasma. The change of graphene 
resistivity upon exposure is confronted with the effect of 
Ar/Il2 treatment in Fig. [3al We see that after Ar expo- 
sure graphene resistivity does not change, remaining in 
the kfi regime and also no D band could be resolved in 
Raman spectra. With these two characterization tech- 
niques we measure no influence of the Ar plasma on the 
graphene devices in spite of strong graphene erosion in 
the case of non-contacted flakes (such flakes are com- 
pletely sputtered after 30 min). This confirms that with 
the chosen plasma conditions no detectable damage is in- 
troduced by Ar ions and that in the flakes with zero self- 
bias the defects detected by Raman spectroscopy come 
only from H adsorbates. These findings also disprove the 
suggestion of Ref. 22] that under exposure to the Ar/H2 
plasma the observed defect band in Raman comes from 
the fragmentation of a water add-layer by high energy 
plasma electrons. If that were the case, we should see 
the Raman band after exposure to Ar in graphene de- 
vices even when the Ar plasma does not introduce defects 
itself. The high energetic plasma electrons from Ar ions 
should similarly fragmentate a water add-layer, which is 
always present in the vicinity of graphene due to the used 
substrate (Si02 is hydrophilic) . Since no Raman band is 
observed after Ar exposure, the Ar plasma does not cause 
graphene erosion and that water layers do not contribute 
to hydrogenation in the plasma process described here. 



E. Relation between the mean free path and defect 
density 
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FIG. 3: (Color online) (a) Resistivity of single (blue dots) 
and double layer graphene (black squares) after se- 
veral exposures to the Ar/H2 plasma. Filled circles 
represent the resistivity at the Dirac point, open cir- 
cles represent the resistivity in a metallic regime (at 
2 X 10^^ cm~^ carrier density). For comparison, filled 
and open diamonds describe the resistivity changes 
in SGL after the Ar plasma exposure. The inset 
presents the exemplary resistivity curve for SLG. (b) 
Mean free path of charge carriers in graphene after 
the exposures. The shaded area indicates the values 
below the length of C-C bond, where the calculations 
of the mean free path loses its physical meaning. 




defect distance Lp [nm] 

FIG. 4: (Color online) Nonlinear correspondence between es- 
timated defect distance and the calculated mean 
free path in single and bilayer graphene. The inset 
presents the measured Id/Ig ratio for SLG and BLG 
devices at three different exposure times. 



Having ascertained that the defects characterized by mean free path to the defect distances L^ extracted 
Raman spectra originate only from H, we can relate the from the Id/Ig ratio. The commonly used Tuinstra- 
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Kroenig experimental dependence [33*] , which relates the 
Id I Iq ratio to the size of graphite nanocrystalinities 
and therefore defect distances, was obtained from X-ray 
diffraction measurements. Estimation of defect concen- 
tration from that relation is inappropriate here, because 
in Tuinstra-Koenig experiment only the edge defects and 
not the whole surface area contribute to the Raman scat- 
tering. We therefore apply a relation established for 
low energy (90 eV) argon ion bombarded graphene from 
Ref. [131, which in the regime measured electronically 
here {Id/Iq <2.5) has a form Id/Ig = (102 ± 2)/L\,. 
The proportionality coefficient was obtained experimen- 
tally for Raman laser wavelength A — 514.5 nm, which 
is close to the one used here (A — 532 nm) and there- 
fore we neglect its possible energy dispersion [s^l- As 
we measure Raman spectra after 3 different exposures 
(see inset in Fig. |4]), the Id/Ig ratios for the exposures 
in between are estimated assuming their linear increase 
in time between the consecutive ratios. The estimated 
defect distance is compared to the electronic mean free 
path extracted from transport measurements in Fig. U) 
We observe a nonlinear relation between the defect dis- 
tances and the mean free path in both SLG and BLG. 
Assuming a parabolic dependence of the mean free path 
from defect distance: I — L^/a we obtain a scattering 
cross section cr of 7 nm for SGL and 4 nm for BLG. 
This confirms that the cross-section for electron scatter- 
ing on the impurity potential is larger than the size of the 
structural disorder caused by this impurity. These scat- 
tering cross-section is roughly the same within the first 
four exposures and then it strongly increases, suggest- 
ing a coalescence of the hydrogenated regions. The lower 
scattering cross-section in BLG supports the theoretical 
predictions that the impurity potential is screened more 
effectively in BLG than in SLG After the last expo- 
sure, the H coverage determined from the defect distance 
Ld is -0.05%. 

As in Ref. |17i] we find that after the Ar/H2 plasma 
exposure the Id/Ig ratio for BLG device is larger than 
that for SLG device (see inset in Fig.H]). This observation 
is in contradiction to the Raman ratios after exposure of 
graphene to atomic H and when other defects are intro- 
duced [13, H^. It is also counterintuitive, as in the bilayer 
the presence of the second graphene layer reduces the rip- 
pling imposed by the amorphous Si02 substrate, which 
should increase the potential barrier for chemisorption of 
H. Also the intensity of the G band in the case of BLG 
should be greater than in SLG, as bilayer resting on a 
substrate can absorb H only on the top layer, leaving the 
layer beneath intact. With the same surface disorder, the 
Id/Ig ratio for BLG is estimated to be 3.5 times smaller 
than for SLG . From that we conclude that the bind- 
ing of H in our process is effectively 4 times larger for 
BLG than for SLG. The observed discrepancy may be 
inherent to the reactivity of ions, the most dominant 
hydrogen-based component in RF plasma, and to their 
dissociation mechanisms at the graphene surface. De- 
tails of this process, together with the exact evolution of 



the Id / Ig ratio with the number of exposed layers, need 
computational verification. Monte Carlo simulations of 
graphite bombarded with H atoms predict that the high- 
est adsorption rate is for H-beam with incident energy of 
5 eV; then in a higher energy range (around 15 eV) the 
H atoms are reflected back from the surface and at even 
higher energies (>30 eV) H atoms are able to penetrate 
through the hexagonal ring and initiate chemical sput- 
tering [s?!- Here the plasma ion kinetic energy ranges 
from 5 eV to 20 eV, which covers both: chemisorption 
and reflection regime for H ions. This may explain the 
somewhat longer exposure times for similar hydrogena- 
tion levels than in Ref. [l3]- This also indicates that the 
efficiency of this process may be still further improved, by 
for example increasing the gas pressure or by increasing 
the RF power. Although the maximum hydrogenation 
limit is not explored here, this plasma technique is ex- 
pected to allow a much higher hydrogen uptake than the 
one reported here (0.05%). 

III. CONCLUSIONS 

In this work we report for the first time the realization 
of graphene hydrogenation in reactive ion etching (RIE) 
system. We study the evolution of the intensity ratio of 
Raman bands Id/Ig and on this basis quantify the in- 
duced disorder. With moderate heating we are able to re- 
verse the hydrogenation to almost initial level, which con- 
firms that the observed disorder in Raman spectra stem 
from adsorbed H. We emphasize here the importance of 
graphene electric potential during the plasma exposure 
to suppress erosion of the flakes. We perform electrical 
studies of single and bilayer graphene after several plasma 
exposures and link them with the amount of the struc- 
tural disorder characterized by Raman spectroscopy. The 
nonlinear correspondence between the mean free path 
and the estimated defect distances is highlighted, from 
which the scattering cross-section for hydrogen defect is 
obtained. We prove that under the chosen plasma condi- 
tions, hydrogenation occurs primarily due to the hydro- 
gen ions and not due to fragmentation of a water add- 
layer by highly accelerated plasma electrons. We also 
demonstrate that by controlling the electric potential of 
the graphene during the plasma exposure we suppress 
the sputtering of carbons in graphene. For that reason 
the hydrogenation level can be precisely controlled and 
reversed. The described hydrogenation process can be 
easily implemented in any RIE system, which we believe 
will stimulate the research of hydrogenated and function- 
alized graphene. 
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